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0.  ABSTRACT 


In  this  paper  some  of  the  problems  of  parametric  estimation  for  samples 
from  distributions,  which  have  decreasing  failure  rate,  are  discussed.  Some 
of  these  distributions  are  obtained  by  mixing  different  exponential  populations, 
others  represent  mechanisms  analogous  to  *  work  hardening*for  each  component  where 
•old  is  better  than  new*.  Sufficient  conditions  are  obtained  that  maximum 
likelihood  estimators  of  appropriately  chosen  shape  and  scale  parameters  exist 
in  both  cases.  The  sample  data  which  are  available,  when  a  decreasing  failure 
rate  distribution  governs  life,  are  usually  censored  with  only  a  few  failures 
observed;  this  limitation  is  dealt  with.  Actual  data,  obtained  from  the 
testing  of  integrated  circuit  electronic  packages,  are  provided  to  illustrate 
these  concepts  and  verify  the  applicability  and  usefulness  of  the  techniques 
described. 
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1 .  INTRODUCTION 


The  exponential  distribution  is  currently  one  of  the  most  widely  used 
models  in  the  study  of  life  lengths,  not  only  because  of  its  early  exposi¬ 
tion,  by  Epstein  and  Sobel  (1953),  and  the  consequent  development  of  the 
associated  tests  for,  and  estimators  of,  its  parameters  but  also  because  of 
its  intuitively  understood  properties,  and  its  theoretical  simplicity. 

The  exponential  family  of  distributions,  serves  as  the  boundary  between 
two  important  classes  of  life  length  distributions,  namely,  monotone 
increasing  or  monotone  decreasing  failure  rates.  These  IFR  and  DFR  families 
respectively,  with  their  distributional  properties,  have  played  central  roles 
in  the  development  of  reliability  theory;  see  the  text  of  Barlow  and  Proschan 
(1975). 

If  each  component  has  a  life  distribution  with  constant  failure  rate, 
the  value  of  which  depends  upon  both  the  variability  in  its  manufacture  and 
the  rigor  of  its  service  environment,  then  life  data  collected  and  pooled  will 
then  exhibit  an  apparent  decreasing  failure  rate.  The  reliability  of  any  such 
component  at  a  time  t  >  0  is  given  by 


where  G  is  the  mixing  distribution  governing  the  failure  rate  variability. 

It  is  known  not  only  that  the  corresponding  hazard  function  -tnR(t) 
will  be  concave  increasing  for  any  mixing  distribution  G,  see  Barlow  and 
Proschan  p.  101,  (1975),  since  the  early  paper  of  Proschan  (1963),  that  there 
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are  instances  when  such  mixing  occurs  naturally  and  therefore  distributions 
with  decreasing  failure  rates  are  of  practical  interest.  Subsequently, 
several  specific  models  of  DFR  distributions  were  discussed  by  Cozzolino 
(1968).  But  parametric  families  of  DFR  distributions  have  received  little 
attention  compared  with  their  IFR  counterparts. 

One  model,  which  is  of  interest  in  its  own  right,  is  based  on  a  gamma 
mixture.  If  for  some  a,  B>  0  we  take  the  density 


G’(A) 


jS-l  e~  VS 

r(a)e“ 


for  A  >  0 , 


it  then  follows  from  equation  (1.1)  that  the  reliability  function  is 

R(t)  »  (l+ter°  for  t  >  0,  (1.2) 

with  cumulative  hazard,  Q  =  -  fcnR , 

Q(t)  =  aJln(l+tS)  ,  for  t  >  0,  (1.3) 

and  failure  rarte  q(t)  =  aB/  (1+tB).  The  density  is  given  by 

f(t)--R'(t)«  - for  t  >  0  .  (1.4) 

(l+tB)a  1 

We  see  by  direct  integration,  or  conditional  expectation,  that 
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E[T]  *  f  R(t)dt 
J  n 


/ 


|“eTa^T) 


1 

I 


2 

B2(a-2)(a-l) 


for  a  >  1 
for  0  <  a  £  1, 

for  a  >  2 
for  0  <  a  £  2, 


and  for  o  >  2 


Var(T)  * 


82(a-l)2(a-2) 


(1.5) 


These  formulae  have  been  given  previously  by  Harris  and  Singpurwalla 
(1968)  and  Barlow  and  Proschan  (1975). 

This  two-parameter  family  seems  to  have  been  introduced  in  this  country 
by  Lomax  (1954)  who  regarded  it  as  a  generalization  of  a  Pareto  distribution, 
which  is  useful  in  business  analysis. 

Kulldorff  and  Vannman  (1973)  and  Vannman  (1976)  have  studied  a  best 
linear  unbiased  estimate  for  the  parameters  of  this  gamma  mixed  exponential 
model.  They  obtained  a  best  linear  unbiased  estimate  of  the  scale  parameter 
B  assuming  that  the  shape  parameter,  a,  was  known  and  in  a  region  restricted 
so  that  both  the  mean  and  the  variance  exist,  namely  a  >  2.  When  this 
restriction  of  a  >  2  cannot  be  met  an  estimate  based  on  a  few  order  sta¬ 
tistics,  which  are  optimally  spaced,  is  claimed  to  be  asymptotically  best 
linear  unbiased  estimate,  and  tables  of  the  weights  as  functions  of  the 
number  of  spacings  are  provided.  In  all  cases,  the  shape  parameter  was 
assumed  known  and  the  sample  was  either  complete  or  censored  at  the  kth  order 
observation.  It  is  contended  that  BLUE  estimates  of  the  shape  parameter  are 
not  attainable. 


Harris  and  Singpurwalla  (1968)  examined  the  method  of  moments  as  an 
estimation  procedure  for  this  same  model  but  again  with  the  shape  parameter 
restricted  to  a>  2,  and  only  for  complete  samples. 

Harris  and  Singpurwalla  (1969)  also  exhibited  the  maximum  likelihood 
equations  for  complete  samples  for  this  model  but  without  any  discussion  of 
computational  methods  for  their  solution, which  was  only  implicitly  defined. 

In  Kulldorff  and  Vannman  (1973)  there  is  a  brief  bibliography  of  results  on 
parametric  estimation  for  the  Lomax  distribution,  with  density  given  by  (1.4), 
under  various  assumptions. 

As  a  consequence  of  the  widespread  adoption  of  integrated  circuitry,  at 
present, life  testing  in  electronic  manufacturing  virtually  always  provides 
incomplete  samples.  This  is  because  tests  are  censored,  because  of  the  expense, 
and  the  paucity  of  failures  due  to  the  high  reliability  of  integrated  circuit 
devices.  Such  service  life  data  cannot  be  adequately  treated  by  any  of  the 
presently  known  statistical  techniques  without  employing  Bayesian  arguments, 
with  their  utilization  of  subjective  information.  This  would  indicate  the 
need  for  an  objective  estimation  procedure,  making  use  of  the  only  type  of  data 
available,  and  without  the  potential  for  bias  inherent  in  Bayesian  priors. 

In  this  paper  the  maximum  likelihood  estimates  are  obtained  for  both  the 
shape  and  scale  parameters,  from  a  class  of  two  parameter  families  (including 
the  Lomax  distribution),  and  necessary  and  sufficient  conditions  are  given  for 
their  existence.  These  estimates  are  derived  for  censored  data  and  a  fortiori 
for  complete  samples,  even  with  a  paucity  of  failure  observations. 
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The  conditions  for  the  existence  of  the  maximum  likelihood  estimates 
apply  even  to  the  case  where  the  mean  and  variance  of  the  life  distributions 
might  not  exist,  e.g.,  in  the  case  of  a  Lomax  distribution  when  0  <  a  <  2  . 
Moreover,  the  estimates  of  the  shape  parameter  a  obtained  from  actual  data 
indicate  it  is  important  to  have  estimation  procedures  valid  in  this  range. 
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2.  DECREASING  FAILURE  RATE  MODELS 

We  now  postulate  the  underlying  physical  processes  which  determine  the 
length  of  life  under  consideration.  Firstly  let  us  suppose  the  quality  of 
construction  determines  a  level  of  resistance  to  stress  which  the  component 
can  tolerate.  The  service  environment  provides  shocks  of  varying  magnitude 
to  the  component  and  failure  takes  place  when,  for  the  first  time,  the 
stress  from  an  environmentally  induced  shock  exceeds  the  strength  of  the 
component. 

If  the  time  between  shocks  of  any  magnitude  is  exponentially  dis¬ 
tributed,  with  a  mean  depending  upon  that  magnitude,  then  the  life  length  of 
each  component  will  be  exponentially  distributed  with  a  failure  rate  which 
is  determined  by  the  quality  of  assembly.  It  follows  that  each  component 
in  service  will  have  a  constant  failure  rate  but  that  the  variation  in 
manufacturing  and  inspection  procedures  causes  some  components  to  be 
extremely  good,  with  a  low  failure  rate,  while  a  few  others  are  bad,  with  a 
high  failure  rate,  and  most  are  in-between. 

Alternatively,  let  us  consider  structures  subjected  to  dynamic 
stresses  of  such  a  nature  that  the  first  stresses  to  which  they  are  sub¬ 
jected  if  not  severe  enough  to  initiate  failure  only  cause  localized 
yielding  and  deformation  thus  affecting  local  stress  relief  and  reinforcement. 
Such  behavior  increases  their  ability  to  withstand  future  stresses  and  could 
be  thought  of  as  "work  hardening".  In  such  cases  older  structures  in  service 
actually  have  greater  resistance  to  fatal  shocks  than  younger  ones;  i.e., 
each  component  in  service  has  a  decreasing  failure  rate.  An  analogous 
behavior  exists  for  increased  strength  or  immunity  in  biological  systems. 
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This  may  be  thought  of  as  "the  older  the  better"  at  least  for  certain 
periods  and  purposes.  Usually  the  failure  mechanism  and  its  interpre¬ 
tation  in  these  cases  is  quite  different  from  that  of  a  mixed  exponential 
model . 

Let  X  be  the  measure  of  the  lack  of  resistance  to  shock  for  a 
component,  the  life  of  which,  say  X^,  will  be  exponential.  If  the  vari¬ 
ability  of  manufacture  determines  the  frequency  of  the  different  X-values, 
which  we  describe  by  a  r.v.,  say  A,  with  distribution  G,  then  T =  is 
the  life  length  of  a  component  selected  at  random  from  those  manufactured; 
it  will  have  a  survival  distribution  which  can  be  written  as  the  condi¬ 
tional  expectation 


EaP[X 


A 


e’XtdG(X)  . 


We  restrict  our  attention  to  mixing  distributions  G  with  unit  mean  the 
choice  of  which  determines  the  distribution  of  life. 

For  the  random  life  length  T  having  reliability  R  we  postulate  a 
cumulative  hazard,  for  some  known  0  <  y  s  1  ,  of  the  form 

r ts 

-S.nR(t)  =  a  I  [q(x)]Y 
J  0 


dx 


where  a  >  0,  0  >  0  are  the  unknown  parameters  for  which,  without  loss  of 
generality,  the  scale  of  X  may  be  chosen  so  that  the  mixed  exponential 
hazard 


q(t)  = 


(2.1) 


/•  00 

is  standardized,  i.e.  q(0)  =  J  X  d  G(X)  =  1 


When  q  and  y  have  been  specified  we  write  T~J(a,B;q,y)  or 
T  ~  J(a,8)  and  we  call  this  the  Afanas'ev  generalization  model  for  life  lengths 
with  unknown  parameters  a, 6. 

It  is  clear  however  that  this  parameterization,  with  y*l,  may  alter  a 
mixed-exponential  distribution  to  one  which  is  not  mixed-exponential  but 
merely  a  decreasing  failure  rate  distribution. 

The  gamma  mixture  of  exponentials  has  been  previously  introduced  in 
equation  (1.2).  Let  us  now  consider  its  Afanas'ev  generalization.  For  a  given 
0  <  y  s  1  we  have  for  its  standardized  hazard  rate  q  , 


q-j (t)  =  (1  +  t)~Y  for  t  >  0  , 


(2.2) 


from  which  we  find  the  corresponding  hazard  function  to  be 

- 

C(1  + 1) 1  “  Y-1 3/(1  -  y)  for  0  <  y  <  1 

<1,(0  *  (2.3) 

$.n(l  +t)  for  y =  1  • 

For  y  =  l>  this  is  a  Lomax  distribution  with  shape  parameter  of  unity. 

It  is  otherwise  a  Weibull  distribution  with  a  location  parameter  of  minus  one 
and  a  shape  parameter  of  less  than  unity.  In  the  Soviet  Union  it  was  introduced 
by  Afanas'ev  (1940)  as  a  distribution  for  fatigue  life. 

We  now  examine  the  distribution  which  results  from  the  mixing  of  two 
exponential  distributions.  It  is  presumed  that  owing  to  occasional  laxity  in 
quality  control  there  is  a  low  probability  p  of  passing  a  component  con¬ 
taining  a  defect  which  can  cause  a  high  failure  rate  X  >  1.  But  there  is  a 
high  probability  of  passing  a  component  having  the  nominal  (low)  failure  rate. 
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which  without  loss  of  generality  we  take  to  be  unity. 

The  reliability  of  this  Bernoul 1 i -mixed  exponential  population  is 


R(x)  =  pe  ^x  +  (l-p)e”x  for  x  >  0  . 

Whence,  we  find  the  standardized  hazard  to  be,  for  x  >  0 

QjU)  =  -£np  +  Xx  -  £n(l  +  re  vx) 

u 

with 

r  =  ^  ,  v  =  — -  >0,  y  =  p(X+r). 

We  shall  write,  respectively, 

T  ~  J^(a,B;y)  or  T  ~  Ogfa.Sip.X) 


(2.4) 


(2.5) 


whenever  the  standardized  hazard  ,  is  defined  in  (2.3)  or  as 
defined  in  (2.4).  Also  without  further  mention  we  shall  use  the  same  sub¬ 
script  to  denote  other  functions  associated  with  these  cases. 

We  now  comment  that  it  is  possible  to  introduce  an  alternative  parameter¬ 
ization  to  that  of  (2.1),  namely 

_  3:  q  ( tg ) 

R(t)  =  e  ®  for  t  >  0;  a,B  >  0,  (2.6) 

where  Q  is  a  known  cumulative  hazard.  The  advantage  is  that  in  this  form 
we  see  clearly  that  B  +  0  implies  R(t)  e  .  Thus  the  limiting  case  is 
exponential  life.  In  either  case  one  sees  that  8  measures  the  degree  of 
departure  from  a  constant  failure  rate  in  terms  of  Q  and  determines  the  rate 
of  improvement  in  expected  life  with  use.  This,  in  a  sense,  is  a  reliability 
growth.  The  disadvantage  is  that  B  in  (2.6)  is  now  no  longer  a  scale 
parameter  and  the  likelihood  equations,  to  be  obtained  subsequently,  seem 


somewhat  more  complicated.  However,  the  transformation  (o,S)  -*•  (a/B,B)  is 
a  1-1  mapping  of  the  positive  quadrant  into  itself  so  that  any  maximum  like¬ 
lihood  estimates  for  one  parameterization  could  be  immediately  transformed  to 
the  other. 
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3.  A  CLASS  OF  TWO  PARAMETER  FAMILIES  OF  MIXED  EXPONENTIALS 

We  shall  postulate  a  class  2  of  concave  hazard  functions  (or  decreasing 
hazard  rate  functions  correspondingly). 

Each  element  of  this  class  will  generate  a  two  parameter  family,  which 
is  a  subset  of  the  DFR  distributions.  Reliability  will  be  of  the  form,  where 
Q  t  2  is  a  given  hazard  function, 

R(t)  =  t  >  0;  a, 8  >  0, 

here  q  =  Q ’  is  a  decreasing  hazard  rate  which  is  twice  differentiable, 
standardized,  i.e.  q(0)  =  1  ,  and  satisfies 

1°  ip  =  xq  is  increasing, 

2°  q  is  log-convex, 

3°  C  =  1  +  xq'/q  is  non-negative,  does  not  exceed  unity 
and  tends  to  a  limit  in  [0,1]  as  x  -*■  °°. 


The  question  arises,  "Where  did  such  assumptions  come  from  and  what  dis¬ 
tributions,  if  any,  satisfy  them?"  One  sees  immediately  that  if  q  €  2  then 
qY  €  2  for  any  y  €  (0,1]  so  that  2  is  closed  under  fractional  powers  of  the 
failure  rate.  The  Afanasev  generalizations  of  the  Lomax  distribution  has 
failure  rate  q-|  given  by  (2.2).  One  sees  that  q,  is  decreasing  and  tnq^ 
is  convex.  Moreover,  one  checks  easily  that  ^  is  increasing  and  in  this 


case 


5,(0  »  [l  +  (l-Y)t]/(l+t) 


(3.2) 


1  w" 
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is  decreasing  with  (®)  =  1-y  .  So  in  this  archtypical  example  assumptions 
1°,  2°,  and  3°  are  met. 

In  the  Bernoulli  mixture  of  the  two  exponentials  we  find  from  (2.4) 


rv 


q*<X)  ‘  5  '  '  (e~v*+r) 


2  — vx 

«2<‘>  ■  7.-11. 


(e  +r) 


and  since  r  >  1  , 


qij(x)  =  v3re~vx(r-e"vx)/  (e'vx+r)3  >  0  , 


(3.3) 


(3.4) 


thus 


C2(x)  =  1 


xrv2e”vx 


( Ae~vx+r) (e"VA+r) 


-vx. 


so  that  we  see  52(°°)  *  ^  * 

We  next  show  that  is  an  increasing  function.  To  see  this  we  note 
that  \p^(x)  >  0  for  all  x  >  0  iff,  after  simplification, 

-ux  2  vx 

Ae  +  r  e  >  vryx-r(l  +A),  for  all  x  >  0. 

VX 

Let  y  =  e  then  the  inequality  above  becomes  clearly  true  for  1  sy<yQ  , 
where  £n  yQ  *  (1  +  A) /(I -A). Let  y  =  yQt  for  t  >  1,  then  the  inequality  to  be 

proved  becomes 

~ ^  +  r2y  >  vry  for  all  t  >  1 . 

But  the  right  hand  side  is  maximized  at  t=e  ,  so  it  is  sufficient  to  have 

yQ  >  vy/re,  which  is  implied  by 

D  < - 1 - —  A  -1 

1  +  (A-l)e2/^'^  X(X  +  1) 


■ 
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We  now  must  show  that  inq  is  convex.  To  prove  this  it  is  sufficient 
to  show  that  qgd £  ?  (qp  •  But  substitution  and  simplification  from  (3.3) 
and  (3.4)  into  the  equation  above  shows  that  this  inequality  is  true  for  all 
x>0  iff  ri  /Y  ,  which  is  equivalent  with  p  :  (1  +  /X)"^  .  We  spare  the 
reader  further  details. 

We  claim  that  these  Inequalities  are  virtually  always  true  in  practice 
since  if,  for  example,  X  should  be  as  high  as  10,  then  the  probability  of 
passing  such  a  bad  component,  with  a  failure  rate  ten  times  the  nominal  design 
rate,  must  not  exceed  .082.  This  would  seem  to  be  a  reasonable  assumption, 
at  least  for  firms  that  intend  to  remain  in  business.  Furthermore  this 
demonstration  shows  that  not  all  decreasing  failure  rate  distributions  satisfy 
our  assumptions. 

Of  course  not  all  DFR  distributions  can  be  obtained  by  an  exponential 
mixture  such  as  given  in  (1.1).  For  an  exact  description  of  the  extreme  points 
of  the  class  of  DFR  distributions,  see  Langberg  et  al  [9].  Moreover  only 
particular  exponential  mixing  distributions  thought  to  be  of  practical  interest 
will  concern  us  here,  along  with  their  corresponding  Afanas'ev  generalizations. 
Thus  we  will  examine  only  a  subclass  of  the  DFR  distributions. 

We  now  ask  what  are  sufficient  conditions  on  the  mixing  distribution  so 
that  assumptions  1°,  2°  and  3°  will  be  satisfied  for  any  Afanas’ev  general¬ 
izations  of  mixed  exponential?  In  the  following  discussion  we  adopt  a  notation 
omitting  the  limits  of  integration  when  they  extend  from  0  to  ®  . 
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Theorem  1 :  If  the  mixing  distribution  6  is  such  that  the  function 
K(x,y)  *  G(x/y)  for  x,y  >  0 


is  totally  positive  of  order  2,  then  1°  is  satisfied. 

Proof:  Making  a  change  of  variable  in  the  definition  we  see 


4>(t)  *  Jye~ydG(y/t)  y/je'y dG(y/t). 


Upon  integrating  numerator  and  denominator  by  parts,  we  obtain 

J G(y/t)ye"ydy 


<l>(t)  +  1  = 


J G(y/t)e"ydy 

Let  tj  >  t2,  we  must  show  t2) -  This  is  true  iff 


0  < 


jG(y/tj  )ye'ydy,  J*G(y/t1  )e-ydy 
/G(y/t2)ye'ydy,  jG(y/t2)e"ydy 


By  applying  the  basic  composition  formula  of  Karlin,  see  e.g.,  p.  100, 
Barlow  and  Proschan  (1975),  to  the  right-hand  side  above,  it  becomes  equal  to 


// 
*1  ‘ 


G( t? / yi ) »  G(tj/y2) 

y^  y1  ,  e  yl 

X 

G(t2/y2) 

y2e  »  e 

^  dy2 


Clearly  the  second  determinant  is  negative  and  that  the  first  is  negative, 
by  definition  of  TP-2,  can  be  seen  by  setting  x^  =  t2,  x2  =  t.| .  QEO. 

The  hypothesis  of  this  theorem  has  a  relation  with  Polya-Frequency 


functions  of  order  2,  if  it  were  expressed  as  the  difference  of  the  logar 
ithns  rather  than  as  a  ratio. 
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We  now  have  the  important  second 

Theorem  2:  For  any  G-mixed  exponential  distribution,  the  failure  rate  q 
is  a  log-convex  and  2°  is  satisfied. 

Proof:  From  equation  (2.1)  we  find 


Unq)’(t)  =  q(  t) 


and  hence 


(£nq)"(t)  ■ 


J*A2e_AtdG(A) 
j Ae"AtdG(A) 

J*A2e"AtdG(  A) 

Ae"AtdG(A) 

Je"AtdG(A) 

*  e"AtdG( A) 

J 

J  A3e~AtdG{A) 

jA2e"AtdG(A) 

I"  Ae"AtdG(  A) 

j Ae'XtdG(A) 

(3.6) 


It  is  sufficient  to  show  that  Unq)"  >  0  .  To  see  the  first  difference 
is  positive  we  fix  t  and  set 

dH(X) 


/' 


e~XtdG(A) 


e‘XtdG(A) 


(3.7) 


H  then  becomes  a  distribution  function  and  the  first  difference  can  be 
written 


jA2dH(A)  -  [  J*AdH(A)]2  , 


which  is  a  variance  and  is  positive.  The  second  difference  is  handled 
similarly.  QED. 


We  now  examine  the  behavior  of  ?  for  a  mixed  exponential,  in  partic¬ 
ular  its  asymptotic  behavior  determined  by  the  behavior  of  the  mixing  distrib¬ 
ution  G  near  0. 
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We  do  not  consider  as  reasonable  for  our  application  mixing  distributions 

x-2 

which  are  of  infinite  order  at  the  origin,  such  as  G(x)  *  e  for  x  >  0  . 

We  admit  for  consideration  only  distributions  G  which  are 
(I)  discrete  in  some  neighborhood  of  zero,  or 
(II)  there  exists  <  >  0  such  that  x-,cG(x)  -*•  c  >  0  as  x  +  0  . 


Theorem  3.  For  any  G-mixed  exponential  distribution,  ?(0)  =  1,  and  c  is 
initially  decreasing.  Since  limit  exists  we  have  c(®)  2[0,1]  and: 

If  G  is  of  type  (i)  then  c(®)  »  1 
If  G  is  of  type  (ii)  then  c(«)  =  0. 


Proof:  Let  us  consider 

,  .  rlt)  =  q'(t)  =  Jx2e'x d6( x/t)  _  fxe~xdG(x/t) 

qU)  jxe‘x  dG(x/t)  j  e"xdG(x/t) 

which  we  see  can  be  regarded  as  a  variance  divided  by  a  mean  and  so  is 


positive  always. 

We  see  c  initially  decreases  linearly  with  slope  q'(0)  since 

C'(0)  =  lim  =  q 1  (0)  . 

x  i  0  x 

We  now  consider  the  limit  as  t  +  <*>.  Integrating  by  parts,  above,  we 


obtain 


*  J(x/t)"K  G(x/t)  x<d(x2e'x) 
J(x/t)’K  G(x/t)  xK  d(xex) 


-  /(x/t)~<G(x/t)  x<d(xe~x) 

J(x/t)‘KG(x/t)xK  de'x 


r(<  +  D 
K<T 


as  t  +  *  . 


The  alternative  when  G  is  discrete  is  verified  by  a  Tauberian  argument 
which  is  similar.  QED. 

Comment:  The  exact  conditions  on  the  mixing  distribution  G  which  insure 
that  1  -  c  is  unimodal  are  not  easily  seen  in  general.  Such  behavior  can  be 
easily  checked  in  the  instances  given. 
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Theorem  4.  For  any  Afanas‘ev  generalization  of  a  mixed  exponential,  say 
q  =  q q  for  some  where  the  standardized  q^  is  given  by 

(2.11)  with  the  mixing  distribution  of  type  (i)  or  (ii)  then  there 
exists  a  v  €  [0,1]  such  that  xv  q(x)  -*■  a  >  0  as  x  +  ». 

Proof: 

Let  us  first  suppose  that  v=l  so  that  q  is  a  mixed  exponential 
hazard  rate.  We  claim  that  v *  0  or  1  suffices  accordingly  as  G  is  of 
type  (i )  or  of  type  (ii) . 

If  (i)  then  for  some  0  <  Aq  <  €  ,  Pq  >  0 

R(t)  =  +  Je’Xt  dG€(A) 

with  the  obvious  interpretation  of  Gf  with  support  inside  (€,»)  ,  in 
which  case 

q(t)  = 

as  t  •+  <*>  .  Thus  v=0  suffices. 

If  (ii)  then  for  some  k  >  0  ,  x'K  G(x)  c  i  0  as  x  0  and 
from  (3.5),  after  integrating  by  parts,  we  have 


Again  multiplying  numerator  and  denominator  x-<  we  have  as  x  -*•  ® 

*(x)  *  Jlx/t)l  G.(t/_x)  tk,d(te~t)  .  K  _ 

W  G(t/x)  tk  de_t 

Thus  v  *  1  suffices. 

We  next  consider  a  failure  rate  which  is  a  fractional  power  of  a  mixed 

exponential  for  which  x  qQ  •*  c*0  for  6  =  0  or  1.  Then  v  *  6y  will 

suffice  since  5 

xVq0Y  *  (x  q0)Y  -  a  >  0  as  x  - 


P0  V  JXe' 


r(x-  xn)t 


dG€(A) 


P0  +  Je‘(X~X 0^  dG6(A) 


X0  *  0 


00 


.  QED. 


-18- 


In  any  life  length  model  one  is  interested  in  the  distribution  resulting 
when  independent  components  having  such  a  distribution  are  put  into  a  series 
system.  From  the  concavity  of  the  hazard  follows  the 

Remark:  If  components  with  independent  life  lengths  T\  ~  J(o.  ,  6^)  i  =  1,...,M 
are  in  series  then  the  life  of  the  system  T  will  be  in  the  same  DFR  class 
and  satisfies  the  stochastic  inequality 

T  *  min  >  J(Ec^  ,  la.  8^/Sa.  ) 
with  equality  when  8^  =  8  . 

Another  property,  which  one  would  ask  of  any  DFR  model,  is  what  improve¬ 
ment  can  be  made  by  "burning-in"  a  component  with  such  a  life?  In  practice 
it  is  often  assumed  that  as  a  result  of  a  burn-in  period,  surviving  compon¬ 
ents  are  exponentially  lived.  In  fact  hurn-in  tests  are  often  required  in 
certain  electronic  component  manufacturing  with  a  statement  of  the  ultimate 
failure  rate  so  obtained. 

Of  course,  not  all  decreasing  failure  rate  distributions  do  become 
constant  after  some  finite  initial  period,  but  that  is  a  model  which  is  often 
thought  to  be  appropriate.  This  indicates  the  importance  of  the  second  model 
introduced  and  of  its  utility  in  a  determination  of  the  economic  value  of  the 
stochastically  extended  life. 

The  residual  life  T^  of  a  component  with  new  life  length  T  and  a 
burn-in  of  duration  x  =  0  is  the  conditional  life  remaining  after  time  x 
given  that  it  is  alive  then;  i .e.  T^  =[T  -x  |  T  >  x ]  . 

The  residual  life  of  any  G-mixed  exponential  is  again  a  mixed  exponential 
but  with  a  different  mixing  distribution.  The  residual  life  will  have  density 
f(y) a  ce-xydG(Y),  where  c  is  the  normalizing  constant. 
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It  is  easily  seen  that  T  for  all  t  >  0  as  long  as  t  is 

within  the  initial  interval  on  which  q  is  strictly  decreasing.  Usually  the 
burn-in  time  t  is  increased  until  any  incremental  decrease  in  the  residual 
failure  rate  is  not  worth  the  increment  cost  of  testing.  This  point  will 
necessarily  be  somewhat  different  for  each  particular  Q  . 

Remark :  A  burn-in  of  t  units  of  time  on  a  component  with  a  new  life 
T  ~  J<| (a,  B)  will  yield  a  residual  life 

Tt  -  J-j[a(l  tB)1  Y,B/(1  +  tB)]. 

If  T  -  ^(a.B'-P.A)  then  a  burn-in  length  t  will  only  alter  the  pro¬ 
portion  of  high  failure  rates,  viz., 

Tt  ~  J2(a,8:p' .A) 

where  the  altered  proportion  is  given 

P'  =  pe"ATB/(pe"AT6+(l-p)e‘T6) . 
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4.  ESTIMATION  OF  PARAMETERS  USING  INCOMPLETE  SAMPLES 

The  samples,  that  are  obtained  when  components  having  a  DFR  life 
distribution  are  tested,  are  virtually  always  incomplete  in  the  sense  that 
testing  is  stopped  before  all  components  have  failed.  A  datum  on  a  com¬ 
ponent  that  "failure  has  not  yet  occurred  after  a  specified  life"  is  called 
an  alive  time  or  a  run-out.  Samples  containing  such  observations  are  censored. 
Samples  in  which  life  tests  are  truncated  at  some  preassigned  ordered  obser¬ 
vation  occur  infrequently,  in  our  experience,  when  electronic  components  are 
tested. 

It  is  assumed  throughout  this  section  that  we  are  given  a  sample 
t  -  (t^,...,t^,...,t  )  where  t|,...,tk  are  ordered  observations  of  times  of 
failure  while  tk+-j , . . . , tn  are  the  ordered  observed  alive- times,  with 
lfksn.  All  observations  are  obtained  from  testing  components  having  a 
«J(a,B)  distribution  with  Q  e  2  and  y  specified  but  with  unknown  parameters 
a  and  B. 

We  now  introduce  notation  for  the  two  empirical  distributions  respectively, 

call  them  F.  and  F  .  We  set 
K  n 

Fk(y)  =  ^  {#  of  t^  <y  for  i  =  l,...,k} 

and  similarly  for  F^,  and  we  make  the  notational  convention,  to  be  used 
subsequently  for  any  function  g,  that 

g(x)  =  J  g(xt)dFn(t),  g(x)=J  g(xt)dFk(t). 

Some  results  will  now  be  given  on  maximum  likelihood  estimation  of  the 
unknown  shape  and  scale  parameters  in  the  case  of  censored  samples  with 
Q  e  D  and  y  e  (0,1]  specified. 
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Remark  5 :  When  the  scale  parameter  6>0  is  known,  there  exists  a  unique 
m.l.e.  of  a,  ‘ay  a,  given  explicitly  by 

a  =  k/nQ(8).  (4.1) 

This  result  is  not  new.  If  8  is  known  and  Q  given,  concave  or  not,  then 
the  values  yi  =  Q(t,.B)  for  i  =  l,...,n  can  be  calculated.  They  are  the  alive 
and  dead  times  from  an  exponential  distribution  with  unknown  failure  rate  a. 

The  total  life  statistic  divided  by  the  number  of  failure  yields  the  usual 
maximum  likelihood  estimate  of  the  mean  life. 

What  is  obviously  true  is  that  not  any  set  of  n  positive  numbers,  say 
_t  =  (t, , . . . ,  t^ . . . .  ,tn)  with  ks  n  designated  as  failure  times  and  the 
remainder  as  alive  times,  can  be  used  to  estimate  both  unknown  parameters  for 
any  Q  e  D  .  In  some  sense  the  sample  must  be  close  to  what  would  be  likely 
from  such  a  distribution.  We  now  have  for  a  specified  Q  e  2  the  result 
Theorem  6 :  For  given  a,;t  there  exists  a  unique  m.l.e.  of  8,  denoted  by 
8,  and  it  is  given  implicitly  as  the  (positive)  root  of  the  equation 

C(x)  -  ^♦(x)-O  (4.2) 

whenever  _  . 

*o/«0>  =  £  <4-3> 

where  ?  is  defined  in  (3.1)  and  xq  is  the  maximum  value  such  that 
X  is  decreasing  n  (0,x)  and  where  6q  =  sup(c(x)/x:  x?  xq>. 

We  remark  that  (4.3)  is  always  true  if  ?  is  decreasing.  One  notes 
from  (3.2)  that  ^  is  decreasing  with  c-j(«)  =  l-y,  but  1S  initially 
convex  decreasing  then  ultimately  concave  increasing  with 

For  completeness  we  state  theorems  6,  7  and  8  together  with  comments  and 
explanation.  The  proofs  are  deferred  until  later. 
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The  failure  to  fulfill  the  condition  (4.3)  for  a  sample  would  indicate 
that  the  presupposed  choice  of  a  or  q  should  be  reexamined.  That  is  to  say, 
either  the  initial  failure  rate  a  is  not  large  enough  to  reflect  accurately 
the  number  of  failures  observed  or  the  induced  function  5  does  not  decrease 
over  a  sufficiently  long  interval. 

Some  simpler  sufficient  conditions  to  check  are  as  follows. 

Theorem  7:  For  given  a,t,  Q  e  D  there  exists  a  unique  m.l.e.  8  of  6 
which  is  given  as  the  solution  of  (3.6),  if  x,  decreases  on  (0,tg)  and 

- K  =  kya"t  tQ  . 

Here  t.,  is  the  maximum  dead  time  and  t  the  maximum  of  tL  and  t  . 
k  m  k  n 

Theorems  6  and  7  say  that  if  a  is  known  then  a  m.l.e.  of  8  exists 
and,  usually,  is  easily  computed.  Otherwise  the  asymptotic  behavior  of 
determined  from  q,  must  be  accounted  for. 

We  now  repeat  this  discussion  when  both  a  and  8  are  unknown.  It  is 
shown  that  both  m.l.e.'s  exist  whenever  the  sample  satisfies  the  condition 

—  T 

2tt  <  t  .  If  the  sample  fails  to  satisfy  this  condition  then  the  model  may 
not  be  appropriate; either  and  a  constant  failure  rate  model  or  a  convex  failure 
rate  model  may  be  indicated  rather  than  a  decreasing  failure  rate  . 

Theorem  8:  For  a  given  sample  _t,  with  Q  e  2  specified  and  a, 6  both 

/v 

unknown,  the  m.l.e.  of  8,  say  8,  exists  as  the  smallest  positive  root  of 
the  equation 

£(x)  -  $(x)  =  0,  for  a  >  0 

where  i  =  ^//Q,  iff  the  sample  satisfies  the  inequality 

7 

2t  t  <  t  . 

When  S  has  been  determined,  the  m.l.e.  of  a,  say  a,  is  then  given  by  an 


(4.4) 
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analogue  of  equation  (4.1),  namely, 

a  =  k/nQ(§).  (4.5) 

We  now  give  the  proofs  of  the  preceding  theorems. 

The  interchange  of  differentiation  and  integration  in  what  follows  is 
justified  easily,  using  the  dominated  convergence  theorem,  since  the  functions 
are  bounded  by  unity. 

Proof  of  Theorem  6:  As  a  result  of  the  testing,  assume  that  we  have  the  vector 
t  =  ( t-j , . . . , t ^ .... tn )  corresponding  to  the  observed  events 

[T.  =  t.]  for  i  =  1 , . . .  ,k  and  [T„.  >  t. ]  for  i  =  k  +  1  ,. . .  ,n. 

II  »  I 

By  definition  the  log-likelihood  L  is  given  by 

,  k  n 

eL  =  n  f(t.)  IT  R(t . ) . 

i=l  1  i=k+l  3 

Substituting  from  (2.1)  and  taking  logarithms  we  find 

l  =  kS.n(a6)  +  2  tnq(t.B)-  al  Q(t-B) 

1  1  1  1 

J*oo  roo 

£nq(Bt)dF.(  t)  -  an  Q(Bt)dF  (£). 

0  K  JQ  n 

Dividing  by  the  constant  k,  we  write 

L(a,B|t)  =  £na  +  £nB  +  Xnq(s)  -  y-  ^(B)  (4.6) 

making  use  of  the  notational  convention  introduced  earlier.  Thus  for  given 
B  >  0,  the  likelihood  L(*|S»t)  is  concave  on  (0,®)  and  the  m.l.e.  of  a 
exists  uniquely  and  can  be  obtained  from  the  equation  L'(a|B,t)  =  0,  the 
explicit  solution  of  which  is  given  in  equation  (4.1). 
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Si  nee  Q  is  concave  increasing  so  is  the  average  Q.  On  the  other 

r\,/ 

hand  tnq  is  convex,  consequently  so  is  £nq.  Thus  for  given  a  >  0  the 
likelihood  L(gja,t)  is  the  sum  of  a  convex  and  concave  function.  It  can  be 
written,  neglecting  constants  depending  upon  a  or  t,  as 

L(8|a,0  3  tni|)(8)  -  ^  Q(£ ' . 

If  we  define  A( 6 )  =  L'(8ja,tJ  we  see  that  we  must  examine  the  roots  of 

A(x)  3  ;(x)  -  ^*(x). 

To  examine  the  behavior  of  this  difference  A,  we  note  since 
C  =  1  +  x(2,nq) ' 


that 
1‘  (x) 


i: 


i 


tc ' ( tx)dFt ( t)  =|  [xt  (£nq)"(xt)  +  t(Hnq)' (xt)]dFk(t). 


By  assumption  (£nq)"  >  0  but  ( Jlnq ) '  <  0.  Now  we  see  by  Theorem  4, 

lim£'  (x)  =  c'(0)t  <  0. 
x+0 

Thus  for  any  fixed  t  with  x  sufficiently  near  zero  the  function  c'(tx) 
will  be  negative  but  for  larger  values  of  x  the  function  may  become  positive. 
Thus  we  conclude  that  the  sample  function  x,  must  initially  be  decreasing 
and  then  could  become  oscillatory,  but  approaches  c(°°),  in  the  unit  interval. 
Note  also  that  ip  is  the  mean  of  increasing  functions  and  is  increasing. 

Recall  that  Xq  is  the  maximum  value  such  that  c  is  decreasing  on 


6q  =  max{£(x)/x 


x  >  xQ}  ; 


(0,x)  and 
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Thus  c(x)  never  crosses  the  ray  5gX  for  x>Xg.  Since  we  have  required 

that  *O/s0)  ^  £ 

there  can  have  been  only  one  value  §  such  that  A(§)  *0  with  §<Xg. 

This  follows  since  ij>(x)/x  is  a  decreasing  function  and  cannot  have  obtained 

the  value  k/an  for  x  <  6Q.  But  uT( x)  cannot  intercept  c(x)  for  x  >  xQ 

since  \j>  lies  above  the  ray  and  $  below  it.  But  because  ^(0)  =  0  <  t(0)  =1 

and  4>(xg)  >  £(xg)  with  \  decreasing  and  increasing,  there  must  have 

been  exactly  one  intersection  prior  to  Xg.  This  proves  theorem  6.  QED. 

Proof  of  theorem  7.  Since  ?  decreases  on  (0,tg),  it  follows  that  c 

must  decrease  on  (0,x)  where  x  =  tg/t^.  It  is  sufficient  to  show 

q7( x)  2  k/an,  since  there  can  have  been  only  one  intersection  previously  and 

none  afterward.  Now  q(tx)  >  K  implies  that 

m 

/CO 

tq(tx)dF  (t)  >  t  K 

0 

since  q(tx)  is  decreasing  as  a  function  of  t  and  qCt^x)  2  q(tmx)  for 
—  —  k 

i  =l,...,n.  Thus  ip(x)  >  xt  K  =  —  .  This  concludes  the  proof  of  theorem 
7.  QED. 

Proof  of  theorem  8.  Consider  the  likelihood- function  L(a,g|t)  defined 

over  the  positive  quadrant  as  given  in  (4.0).  All  the  stationary  points,  which 

are  determined  by  t_,  can  be  found  from  the  simultaneous  solution  of 
31  31 

—  *  0,  *  0.  This  yields  the  two  equations  in  a  and  6; 

od  <jOl 

t(6)*f-?(6);  a  5(B). 

Combining  these  into  a  single  equation  in  the  unknown  6,  we  are  led  to  seek 
&  as  a  zero  of  the  difference  s(x)  -  4>(x)  for  x  >  0  where  $  = 
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Consider  the  likelihood  as  given  in  (4.6)  as  a  function  of  a, 8  into 
which  we  substitute  the  relation  displayed  above  to  obtain  the  likelihood 
maximized  with  respect  to  o  for  any  value  of  8.  Neglecting  constants 
depending  upon  t,  this  function  can  be  written,  now  using  argument  x,  as 

L(x)  =  7nq(x)  -  £n[()(x)/x]. 

Interchanging  the  order  of  integration  in  Q(x)/x  we  obtain 

Q(x)/x  «  t  j  q(tx)  dW  { t)  =  t  •  q*(x) 

JO  n 


where  Wn  is  the  distribution,  with  density  given  by 


w;(t)  =  [1  -Fn(t)]/t 


for  t  >  0. 


Again  neglecting  constants,  we  have 


2.nq(tx)  dF.  (t)  -  S£n[  q(tx)  dW  (t) 
0  *  JO  " 


=  £nq(x)  -  Hnq*(x)  . 


(4.8) 


Using  a  Maclaurin  expansion  of  q  we  see 


L(x)  =  x  q * ( 0 )  [t  -  t*]  +0(x2) 


as  x  +  0, 


where 


t  =  r t  dFk ( t)  and  t*  =  f  t  dWp(t)  =  ~?/2t 


Thus  L  is  positive  in  a  neighborhood  of  zero  iff  t  <  t  ;  and  since 

xL‘  =  E  -  4  this  difference  will  be  positive  in  a  neighborhood  of  zero  under 


the  same  condition. 
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We  now  examine  the  behavior  of  L  under  the  more  strict  condition 
that  the  sample  distributions  and  Wn  as  defined  satisfy  Fk>Wn- 
It  then  follows  that  q  5  q*  and  since  £n  is  concave,  that  £nq*>£nq. 
Hence  from  (4.8)  we  see  that 

£nq  -  £nq  5  L  5  £nq  -  £nq  . 

Now  one  notices  the  lower  bound  is  always  negative  while  the  upper  bound  is 
always  positive.  But  near  the  origin  L  is  approximately  equal  to  the 
upper  bound  while  for  larger  values  of  its  argument  L  is  approximately 
equal  to  the  lower  bound  and  for  intermediate  values  L  makes  a  smooth 
transition  between. 
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5.  THE  COMPUTATION  OF  §  FOR  ARBITRARY  SAMPLES. 

As  a  matter  of  practical  calculation  we  are  concerned  with  the  smallest 
root  of  the  difference  l  -  4,  in  the  case  when  c  Is  monotone  decreasing. 

Let  us  consider  the  composite  function  f(x)  =  c"^[$(x)]  in  a  neighborhood 
of  zero  with  the  location  of  the  smallest  crossing,  if  there  is  more  than  one, 
of  the  45°  line.  These  are  equivalent. 

1  —  J. 

An  alternative  expression  for  $  is  4  =  xQ  /Q  ,  with  Q  a  convex 

—  _  J. 

function  decreasing  between  t  and  t  •  q(®)  ,  while  "the  smoothed  Q  "  , 

viz.,  Q/x,  decreases  between  the  same  limits  at  a  slower  rate.  It  follows 

that  Q/x  >  Q  >0  so  that  0  s  <J>  5  1 .  The  S>  begins  at  unity,  initially 

decreases  at  a  decelerating  rate  and  tends  ultimately  to  s(«>) .  To  see  this 

note 

♦  <-)  -  11m?'  4  •  11.  ft'-v(xt)',q(txk(£x)OFn(t)  ^ 

00  X-«0  f  j  ' 

J  t'"v(xt)Vq(tx)dFn(t) 

since  by  theorem  4  there  exists  a  v  e  [0,1]  such  that  xvq(x)  -*•  a  *  0,  as  x -►<*>. 
The  composite  function  f  *  z,  behaves  in  a  neighborhood  of  zero  as  a 
contractive  map,  being  initially  greater  than  x,  then  crossing  at  8  are 
one  being  below  x  for  a  range,  perhaps,  thereafter. 

Thus  we  know  that  successive  iterates 

6i+l  =  f(^i )  for  i  =  0,1,2  ... 

will  converge  to  B  >  0  as  long  as  Bq  <  B*»  the  next  larger  zero  of 
£  -4,  if  one  exists.  Otherwise  the  iteration  will  converge  to  zero. 

Moreover  for  the  special  cases  in  (3.2)  the  inverse  t"1  can  be  easily 


found. 
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Our  computational  procedure  based  on  the  sample  is  as  follows: 
Algorithm:  Given  as  failure  times  and  tk+1,...,tn  as  censored 

times  from  a  mixed  exponential  distribution  with  prescribed  Q  e  D  proceed 
as  follows: 

(1)  Compute  the  sample  moments  t,  t,  t  . 

(ii)  If  t  <  2t  •  t,  assume  the  observations  are  from  an  exponential 
distribution  with  failure  rate  X  and  estimate  it  by 

^  lc 

X  =  — 
nt 

~2  —  - 

(iii)  If  t  >  2t  t,  assume  the  observations  are  from  a  Mixed  Exponential 
distribution,  with  prescribed  Q  e  2  . 

Using  the  sample  functions  explicitly  given  by 

n  n  _  k 

4>(x)  =  X  I  t.q(xt.)/  I  Q(xt. )  l(x)  *  1 +x  E[t.q'(xt.)/q(xt.)], 

1  1  1  1  1  1  1  1  1 

we  guess  6g  ,  then  iterate  using  the  inductive  step; 

given  B-,  compute  $(6^)  and  calculate  Bi+1  such 

that 

we  find  B=lim  0. ,  and  compute 

a  -  k/  S  Q(t.B). 

1  1 

The  nature  of  the  intersection  guarantees  that  within  the  region  when 
♦  and  ?  both  decrease  the  iteration  will  rapidly  converge,  with  a  rea¬ 
sonable  first  guess.  When  the  functions  q  and  Q  are  simple  a  small  pro¬ 
grammable  electronic  calculator,  such  as  the  HP-67,  can  be  used  to  obtain 
these  estimates. 


We  now  turn  to  an  examination  of  the  conditions  on  a  sample  which  will 


guarantee  that  it  can  satisfy  (4.4). 

Since,  remarkably  enough,  condition  (4.4)  is  the  same  for  any  Q  z  2 
we  shall  say  a  sample  from  a  DFR  distribution,  containing  k  dead  times  and 
n-k  alive  times,  is  ample  whenever  it  satisfies  (4.4);  i.e.,  the  following 
inequality  is  met, 

k  n  no 

2  Z  t.  E  t.  <  k  it/.  (5.1) 

1  1  1  J  1  1 

One  can  check  that  a  complete  sample  of  failure  times,  i.e.,  with  k =  n, 
will  satisfy  (5.1)  iff  the  sample  standard  deviation  exceeds  the  sample  mean. 

Of  course,  we  know  that  for  any  DFR  distribution  the  standard  deviation  does 
exceed  the  mean,  when  they  both  exist.  Thus  for  most  complete  samples  from 
Mixed  Exponential  distributions,  (which  are  sufficiently  large),  the  require¬ 
ment  (5.1)  should  be  satisfied. 

Remark:  A  truncated  sample  with  k  <  n  failure  times  and  the  remaining 
observations  truncated  at  tQ  will  satisfy  (5.1)  if 

*0  =■'<’*  /hH-*')  *  ‘t-TT-  for  n  sra11' 

To  see  this  note  It..  2  >  (n-k)tg^  and  so  (5.1)  must  hold  if 

t  ^  >  2t  t  +  — -k—  (t)^ 

By  the  quadratic  formula  this  is  equivalent  to  the  inequality  stated.  The 
second  approximation  follows  from  the  first  two  terms  of  the  binomial  expansion. 

Some  studies  of  the  frequency  with  which  samples  from  mixed  exponentials 
satisfy  condition  (5.1)  have  been  made  by  Sunjata  in  an  unpublished  thesis  (1974). 
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6.  THE  CONFRONTATION  WITH  DATA 

We  now  present  some  data  sets  from  two  different  lots  of  flight  control 
electronic  packages.  Each  package  has  recorded,  in  minutes,  either  a 
failure  time  or  an  alive  time.  An  alive  time  is  the  time  the  life  test  was 
terminated  with  the  package  still  functioning.  j 

First  Data  Set 
Failure  times: 

Alive  times: 

Second  Data  Set 
Failure  times:  37,  53 

Alive  times:  60,  64,  66,  70,  72,  96,  123. 

One  checks  that  both  data  sets  satisfy  condition  (5.1)  so  that  both 
parameters  can  be  estimated  in  a  Mixed  Exponential  model.  We  first  compute 
estimates  of  the  failure  rate  using  the  total  life  statistic: 


1,  8,  10 

59,  72,  76,  113,  117,  124,  145,  149,  153,  182,  320. 


i 
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time  t 

in  min 

Data  Set  1 

Data  Set  2 

X  »  .00197 

estimate  of  reliability 

A 

at  time  t,  R1  ( t) 

X  -  .00312 

estimate  of  reliability 

at  time  t,  R2(t) 

6 

.988 

.981 

10 

.980 

.969 

30 

.943 

.911 

50 

.906 

.856 

100 

.821 

.732 

130 

.774 

.667 

If  we 

assume  the  data  are  observations 

from  the  mixed  exponential 

distribution  of  equation  (1.2)  then  using  the  estimation  techniques  derived 
previously  in  this  paper  we  have  the  following  estimates: 
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time  t 
in  min 

6 

10 

30 

50 

100 

130 


Data  Set  1 
a  =  .0453  8  =  1 .03 
estimate  of  reliability 

A 

at  time  t,  R^(t) 

.915 

.896 

.855 

.836 

.810 

.801 


Data  Set  2 
a  =  .420  8  =  .01 

estimate  of  reliability 

A 

at  time  t,  R2(t) 

.976 

.961 

.896 

.843 

.747 

.705 


R(t)  >e" 


.047t 


R(t)  >  e 


- .0042t 


Note  that  these  estimates  are  more  consistent  with  what  the  data  show,  that  is, 
for  at  least  the  first  50  minutes  we  expect  the  reliability  estimate  for  the 
second  set  of  data  to  be  higher  than  the  reliability  estimate  for  the  first  set 
of  data.  Beyond  this  time,  however,  say  at  100  minutes,  the  data  indicate  that 
the  reliability  estimate  from  the  first  set  of  data  should  be  higher  than  the 
reliability  estimate  from  the  second  set  of  data.  Using  mixed  exponential 
estimates  this  is  the  case. 

A  statistical  test  to  determine  whether  the  data  require  a  constant  or 
decreasing  failure  rate  was  run  on  the  data  from  sets  1  and  2.  For  data  set 
1  we  reject  constant  failure  rate  (in  favor  of  decreasing  failure  rate)  at 
the  .10  level.  For  data  set  2  we  cannot  reject  the  constant  failure  rate 
assumption.  In  this  case,  however,  the  constant  failure  rate  estimates  for 
reliability  and  the  mixed  exponential  estimates  for  reliability  are  close. 


7.  CONCLUSION 


If  a  component  has  a  life  distribution  with  an  increasing  failure  rate, 
the  information  necessary  to  estimate  its  parameters  must  contain  failure 
times.  In  practice  this  means  that  if  there  are  few  observed  failures,  within 
a  fleet  of  operational  components,  there  is  little  information  with  which  to 
assess  their  reliability.  If  a  component  has  a  constant  failure  rate  then 
both  failure  times  and  alive  times  contribute  equally  to  the  estimation  of  its 
reliability.  This  study  suggests  that  if  a  component  has  a  life  distribution 
with  decreasing  failure  rate  it  is  the  alive  times  within  the  data  which  con¬ 
tribute  principally  to  the  estimation  of  the  parameters  (and  thereby  to  the 
determination  of  the  reliability)  since  only  one  failure  observation  is  required 
even  to  estimate  two  parameters,  presuming  the  data  is  ample. 

The  usual  justification  for  using  maximum  likelihood  estimates  is  due  to 
their  asymptotically  optimal  properties,  and  to  their  asymptotic  normality. 

The  problem  of  obtaining  the  usual  sampling  distributions  of  the  MLE's  of  the 
parameters  obtained  for  Mixed  Exponential  model  seems  to  be  difficult,  not 
only  because  the  estimates  are  only  implicitly  defined,  but  also  because  samples 
are  virtually  always  censored.  Furthermore,  the  usual  proofs  for  the  asymptotic 
optimality  of  the  MLE's  may  not  apply  when  censoring  is  of  a  general  type  and 
when  only  sparse  failure  data  is  available.  A  useful  asymptotic  theory  must 
be  developed  for  censored  data  sets  of  which  the  life  of  electronic  packages 
of  integrated  circuits  are  an  illustration. 
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